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Abstract. We present high angular resolution interferometric observations of the 3 and 1.3 mm continuum emis- 
sion, and HCO'''(1-0) and SiO(2-l) v=0 lines, obtained with the Owens Valley Radio Observatory millimeter-wave 
array, toward two hot cores (HCs) associated with two well known ultracompact (UC) Hll regions: G29.96— 0.02 
and G31. 41-1-0. 31. These HCs are believed to host young forming massive stars which have been suggested to be 
surrounded by massive rotating accretion disks. The aim of these new observations is to study the structure and 
kinematics of the molecular clumps surrounding the HCs and nearby UC H ll regions at moderately high angular 
resolution. Our observations reveal that the clumps within which the HCs and UC Hll regions are embedded have 
a complex kinematical structure. The total mass of the clumps is estimated to be in the range 1000-3000 Mq, 
consistent with previous findings. Our observations also show compelling evidence that the clump in G29.96— 0.02 
is contracting onto the HC position, suggesting that the accretion process onto the massive young stellar object 
embedded in the HC is still ongoing. In these objects the kinematical structure that we observe is also compatible 
with the presence of a massive rotating disk within the HC, even though we cannot prove this suggestion with our 
data. The case of G31. 41-1-0. 31 is more complicated, and our data, although consistent with the presence of an 
inner disk and an infalling envelope around it, do not have the required spatial resolution to resolve the different 
structures. 
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1. Introduction 

The importance of the study of high-mass star forma- 
tion derives from the close interaction between such stars 
and the interstellar medium: they affect the ISM evolu- 
tion and, consequently, the galactic evolution in general. 
In spite of this, until the last few years the study of star 
formation was mostly concentrated on low-mass stars, be- 
cause of the difficulties in investigating young massive 
stars, less numerous, more distant, more embedded and 
with much shorter evolutionary timescales. The search for 
the sites where the first phases of high-mass star forma- 
tion occurs has been increasing recently, because of the 
advent of aperture synthesis mm-interferometers (Kurtz 
et al. pOOOl). 



Due to the large luminosity and Lyman photon emission, 
a newly born massive star heats the dust and ionises the 
gas in the surrounding envelope and can be detected in 
the far-infrared (FIR) as a luminous IRAS source and in 
the centimetric radio continuum as an ultracompact (UC) 
Hii region (Wood & Churchwell |l989a{ |l989"b| ). However, 
these are signposts of already-formed massive stars; in or- 
der to study the more elusive evolutionary phase of main 
accretion one has to devise suitable target selection proce- 
dures. Several different approaches have been explored to 
obtain lists of candidate massive accreting protostars (see 
the review by Kurtz et al. 2000): high- luminosity IRAS 
sources with peculiar colours (e.g. Molinari et al. 1996; 
^99§), H2O masers without associated radio continuum 



High-mass stars ignite hydrogen burning and reach 
the main sequence while still accreting matter (Palla & 
Stabler 1993| ). Thus young massive stars are still deeply 
embedded in a dense gas and dust envelope which impairs 
the detection at optical and near infrared wavelengths. 



emission (Tofani et al. |1995| ; Cesaroni et al. |1999bD , com- 
pact molecular cores close to UC Hii regions (Cesaroni 
et al. 1991; 1994a). The latter approach is supported by 
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the observational result that massive stars are usually 
found in rich clusters and associations: it is thus likely 
to find a forming massive star l ying c lose to a young 
one. In this way Cesaroni et al. (1994a) identified com- 
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pact molecular cores close to, but generally not coincident 
with, the target UC Hii regions. These are characterized 
by bolometric luminosities exceeding 10^ Lq, small di- 
ameters (< 0.1 pc) and kinetic temperatures over 50 K, 
and were then called hot cores (HCs). Due to their high 
luminosity and temperature, HCs have been suggested 
to contain very young massive (proto-)stars. In fact, a 
few of them contain UC Hii regions (e.g. GlO.47-1-0.03, 
see Cesaroni et al. 199^ ), which witness the existence 
of embedded early type stars. Moreover, two well known 



HCs, IRAS 20126+4104 (Cesar oni et al. p97| ; [l999a) and 
G9.62+0.19-F (Hofner et al. |l996b| ; Testi et al. |l998; 



2000|) , show evidence of an embedded massive young stel- 
lar object in a phase preceding the formation of a classical 
UC Hii. 

Two of the best stu died classical HCs from the list of 
Cesaroni et al. ( 1994a ) are those close to the UC Hii re- 
gions G29. 96-0.02 and G31.41+0.31. In the former region 
the molecular core lies in front of the cometary UC H ii 
region, approximately 2'.'6 west of the 1.3 cm continuum 
peak. In G31.41+0.31, faint NH3 emission is detected to- 
ward the extended halo around the UC H 11 region, while 
the NH3 HC peak is ^A" to the southwest of the ccntimet- 
ric radio continuum peak. High angular resolution (^0.4") 
Very Large Array (VLA) observations of the NH3(4,4) in- 
version transition (Cesaroni et al. 1998) showed resolved 
emission from both HCs, with a deconvolved size of ^1'.'4 
in G29. 96-0.02 and -I'.'S in G31.41-I-0.31. Both cores ap- 
pear shghtly elongated: G29.96— 0.02 in the E-W direc- 
tion and G31.41+0.31 along a SW-NE axis. A velocity 
gradient of approximately 2-3 km arcsec"^ over 3" is 
found along such directions. Similar correlations are found 
in the CH3CN(6-5) transition for both G29. 96-0.02 and 
G31. 41-^0.31, although the velocity gradients in CH3CN 
(~15 km arcsec" ^) are larger than those observed in 
NH3 (Cesaroni et al. 1994b ; Hofner et al. in prep.). Also, 
in both sources the derived kinetic temperature shows a 
peak toward the center. All these features have been inter- 
preted as evidence of disk-like structures rotating about a 
hot, massive central object: the (proto-)star. 

The presence of a deeply embedded massive 
(proto-)star is also indirectly suggested by the presence of 
several other signposts, such as H2O and H2CO, CH3OH 
masers toward the hot ammonia core peak of G29.96— 0.02 
(H ofner fc Churchwell |l9964 Pratap et al |l994t Walsh et 
al. |1998D and H2O masers toward G31.41+0.31 (Gaume 
& Mutel |1987]). 



Recently, Pratap et al. ( |1999[ ) investigated the molec- 
ular gas surrounding G29.96— 0.02 using the BIMA mil- 
limeter wave interferometer with moderate angular reso- 
lution (~10"). They found ^^CO and high-density tracer 
(CS and CH3CN) emission; while the former molecule is 
found in an extended cloud surrounding both the HC 
and the UC Hll, the latter peak toward the HC posi- 
tion. Hence, CH3CN and NH3 emission trace the inner- 
most and warmer molecular gas within the HC itself, while 
^•^CO traces the large-scale molecular cloud in which the 
HC, the UC Hii region and the surrounding near infrared 



young stellar cluster are embedded. Pratap et al. (1999) 
also detected a compact core of CH3OH emission, ~5" 
south west of the HC. This second core does not show the 
typical features of HCs (compact NH3 and CH3CN emis- 
sion) and Pratap et al. suggest that it may be associated 
with a molecular outflow from the HC exciting source or 
with a low-mass star forming outside it. 

In order to study the relationship between the HC and 
the UC Hii region, and the kinematics of the gaseous 
clump in which these are embedded, we decided to obtain 
high angular resolution observations of the HCO'''(1-0) 
and SiO(2— 1) v=0 lines, as already done in similar regions 
by other authors (Acord, Walmsley & Churchwell 1997; 
Cesaroni et al. 1997] ; Wilner et al. 1996). These tracers 
are expected to provide information on the innermost part 
of the cloud and the kinematics of the molecular gas at 
the interface between the HC and the UC Hii region. 
The main goal of the observations was thus to relate the 
observed velocity field within the HC to the kinematics 
of the molecular gas in the inner portions of the cloud. 
Additionally, the SiO(2— 1) v=0 transition, which can be 
enhanced in shocked regions, was exp ected to clarify the 
CH3OH observations of Pratap et al. ( 1999 ). 

Throughout this paper we will assume a distance of 
6 kpc for G29 .96— 0.02, based on the recent estimate of 
Pratap et al. ( |l999D , and 7.9 kpc for G31.41-f0. 31, b ased 
on the kinematic estimate of Churchwell et al. ( 1990 ). 

We present our new millimeter interferometric observa- 
tions and data reduction in Sect. |^. In Sect. |^ the observa- 
tional results are presented for the millimeter continuum 
and molecular lines separately. We discuss the implica- 
tions of our results on the structure of the two regions 
and the nature of HCs in Sect. ^. In Sect. || we summarize 
the main results of our study. 



2. Observations and data reduction 

High resolution interferometric observations were per- 
formed in the period October-December 1997 using the 
Owens Valley Radio Observatory (OVRO) millimeter- 
wave array, located near Big Pine, California. The array 
comprises six 10.4 m dishes that were employed in three 
different configurations yielding baselines in the range 15- 
220 m. The phase centers used were: 01950 = 18^43™2??1, 
,^1950 = -02°42'36'.'4 for G29.96-0.02 and aiggo = 
18h44™59!l, (S1950 = -01°16'07'.'3 for G31.41-h0.31. The 
half-power beamwidth of the OVRO antennas is ~84" at 
88 GHz, and ~32" at 233 GHz. Each antenna is equipped 
with cryogenically cooled SIS receivers that offered typical 
system temperatures of ^300 K at 88 GHz, and ^1200 K 
at 233 GHz. The 3.4 mm and 1.3 mm continuum obser- 
vations employed an analog correlator for a total band- 
width, at each of the two frequencies, of ~2 GHz. The 
flexible digital correlator was configured to observe simul- 
taneously the HCO+(1-0) and SiO(2-l) v=0 transitions 
at 89.2 and 86.8 GHz respectively. For both molecules we 
used a ^105 km/s wide band with '^1.7 km/s resolution. 
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The spectrometer modules were centered at '^98 km/s, 
close to the systemic velocity of both sources. 

Complex gain calibration was ensured by frequent ob- 
servations of 1741—038. The typical observing cycle alter- 
nated scans on G29. 96-0.02, G31.41-h0.31 and 1741-038 
with a period of ~30 min. 3C273 was used as passband 
calibrator and Uranus and/or Neptune were used to set 
the flux scale. The expected flux calibration uncertainty 
is '^20%. Data calibration and editing has been perfor med 
using the MMA software package (Scoville et al. 1993 ). At 
230 GHz data could be salvaged only for a small subset 
of the observations, due to the poor phase stability and 
relatively high water vapour content of the atmosphere 
during most of the observations. This resulted in a poor 
{u,v) coverage of the final 233 GHz dataset, a final noise 
higher than normally expected and a much reduced sen- 
sitivity to the extended structures. The calibrated {u,v) 
datasets were exported to the NRAO-AIPS package for 
imaging and analysis. All the maps presented in this paper 
were obtained with the AIPS IMAGR task using natural 
weighting of the (u, v) data. 

In the continuum maps, the final synthesised half 

power beam width (HPBW) is 5'.'9 x 3'.'7, with P.A. 8, 

at 88 GHz, and 4'.'9 x 2'.'7, with P.A. 3.5, at 233 GHz, 

for both sources. The rms level is 10 mJy/beam at 88 GHz 
for both sources, 40 mJy/beam for G29.96— 0.02 and 160 
mJy/beam for G31.41-f0.31 at 233 GHz. In the line cubes 
we obtained 5"9 x 3'.' 7 HPBW and 50 mJy/beam rms level 
for all lines. The P.A. is -10.5 for the HCO+ and -16.1 for 
the SiO for G29.96-0.02; -9.2 for the HCO+ and -14.9 
for the SiOline for G31.4H-0.31. The conversion factor 
from flux density to brightness temperature in the syn- 
thesised beam is ^10.4 K (Jy/beam)~^ at 88 GHz and 
-2.43 K (Jy/beam)-i at 233 GHz. 

2.1. Continuum maps 



Our millimeter continuum maps are sensitive to both the 
free-free emission from the ionised gas inside the UC 
Hii regions and the thermal dust emission. The latter 
is mainly confined to the HC region (Olmi et al 



G29.96-0.02 



1996 



Patrap et al. 199£ ). In order to disentangle the HC con- 
tinuum emission from that of the UC H ii region, we need 
to estimate the contribution of the free-free emission at 
millimeter wavelengths and subtract it from our contin- 
uum maps. To this purpose we have used the 24 GHz 
VLA observations of Cesaroni et al. ( |1994a| ; |1998| ). We 
re-imaged the calibrated (u, v) data using the same (u, v) 
range and restoring beam as the 88 and 233 GHz OVRO 
observations. These maps were then scaled assuming opti- 
cally thin free- free emission {Fi, oc v~'^-^). In Figs. ^andU, 
we show for each source and each frequency the millimeter 
continuum maps as observed at OVRO (top panels), the 
free-free emission maps obtained by extrapolation of the 
VLA data (central panels), and the residual "dust" emis- 
sion maps after removal of the free-free emission from the 
millimeter observations (bottom panels). 
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Fig. 1. Top panels: 88 and 233 GHz millimeter continuum 
maps observed at OVRO towards G29.96— 0.02. Central 
panels: expected free-free contribution computed from the 
24 GHz VLA data of Cesaroni et al. ( |l994at |1998D us- 
ing the same (u, v) range and restoring beam as in the 
top panels. Bottom panels: residual "dust emission" maps 
after subtracting the central panels from the top panels 
maps. The filled triangle marks the ammonia HC peak 
position. The beamsize at each frequency is shown on the 
top panels. 



Figures |l| and || show that the procedure described 
above allowed us to compute a fairly accurate estimated 
map of the free-free emission contribution at millimeter 
wavelengths. The subtraction of the extended UC H ii is 
successful in all cases except for the 233 GHz observations 
of G29.96— 0.02. In this case, the poor (m, v) coverage of 
the millimeter map does not allow a proper reconstruc- 
tion of the extended emission, and even if the free-free 
map was calculated using the same (u, v) range as the mil- 
limeter map, the much finer VLA sampling of the {u, v) 
plane allowed us to properly map part of the extended 
emission. When the subtraction of the calculated map is 
performed, the free- free emission from the UC H ii is over- 
subtracted, producing a large negative hole at the position 
of the UC H II. On the other hand, also in this case, the 
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Table 1. SED fit parameters 
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Fig. 2. As in Fig. for G31.41+0.31. 



position of the excess "dust" emission at 233 GHz is co- 
incident with the NH3 core and the 88 GHz "dust" peak. 
Toward G31.4f-|-0.31 the free-free contribution to the mil- 
hmeter emission is lower: at 88 GHz, the observed contin- 
uum peak is twice the value calculated from the 1.3 cm 
map, whereas at 233 GHz the maximum of the expected 
free-free emission is less than the 1 a noise level in the 
observed continuum map. In this case, the observed flux 
at 233 GHz is essentially thermal dust emission. However, 
as for G29. 96-0.02, the total continuum flux at 233 GHz 
measured in our interferometric maps is less (a factor of 
^ 3 for G31. 41-1-0. 31) than what expected from previous 
single dish observations (see Fig. ^), probably due to the 
spatial filtering of extended emission (see also Sect. 3.1). 



2.2. Line maps 

Continuum subtraction from the line data cubes was per- 
formed both in the (m, v) plane using the AIPS task 
UVLIN, selecting line free channels prior to image de- 
convolution, and in the image plane after cleaning, by 
subtracting the appropriate continuum images. The two 
techniques offered comparable results. The analysis pre- 



source 


d 


run 




logiVLy 


Tdust 


/3 




(Kpc) 


(pc) 


(K) 


(ph./s) 


(K) 




G29.96-0.02 


6 


0.5 


6500 


48.90 


30 


2 


G31.4H-0.31 


7.9 


0.5 


8000 


48.75 


37 


1 



sented in this paper uses the continuum subtracted cubes 
produced with the second method. 

3. Results 

3.1. Continuum 

In Fig. H we show the spectral energy distribution (SED) 
for G29. 96-0.02 (top) and G31.41-H0.31 (bottom). The 
flux density measurements that we used are a compilation 
of single dish and interferometer observations (see refer- 
ences m Fig. 1). The data are heterogeneous in both angu- 
lar resolution and (it, v) coverage. The SED is that typical 
of UC Hii regions, see also Wood & Churchwell ( |l989bD , 
dominated by free-free emission at centimeter wavelengths 
and by thermal dust emission at millimeter and far- 
infrared wavelengths. 

To give an estimate of global parameters, such as the 
supply rate of Lyman photons required to mantain the 
UC H II region ionised and the average temperature of the 
dust, we used a simplified two component model: a spheri- 
cal and homogeneous H 11 region, plus a single temperature 
grey-body dust component. This simplified model is able 
to accurately reproduce the radio to far infrared portion of 
the SEDs shown in Fig. |[ In Table |l| we report the values 
of the fixed parameters of the models: distance [d) and 
electron temperature (Thii) of the UC Hii regions. For 
G29.96— 0.02, we used the values estimated by Pratap et 



al.( 199£ ); for G31. 41-1- 0.31, the distance is that estimated 
by Churchwell et al. ( 1990 ), whereas the temperature is 
a mean value for a typical Hii region. Also shown in the 
table are the values of the parameters derived from our 
model fits: UC Hii radius (rfjii), Lyman photons supply 
rate (TVLy), dust temperature (T^), and dust opacity index 

From the SEDs reported in Fig. |3[ we can clearly see 
that our 233 GHz observations are well below both the sin- 
gle dish points and the model fits, expecially in the case of 
G29.96-0.02, where the contribution of the UC Hii free- 
free emission at such high frequency is still a consistent 
fraction of the total emission. This is mainly due to the 
sparse (u, v) coverage of our synthesis at this wavelength, 
and the inability to properly image extended structures. 
Hence, as discussed in Sect. 2.1, our 233 GHz observations 



cannot be used to infer the properties of the large scale 
structures in these regions, but offer only an estimate of 
the high frequency emission of the HCs themselves. 

The values of Tdust and (3 are average values for the 
dust in the envelope surrounding the HCs, not the appro- 
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Fig. 3. Spectral energy distribution for G29.96— 0.02 and G31.41+0.31. The solid lines are fits of the SED assuming 
a simplified two component model: free-free emission from a spherical, homogeneous and dust-free Hli region (dotted 
line), and grey-body thermal dust emission (dashed line). On the right we indicate the reference for each of the 
measurement reported in the plots. 

Table 2. Observed parameters of the continuum (dust and free- free) emission 
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priate values for the HCs since these will only contribute 
to a small fraction of the single dish submillimeter and far 
infrared fluxes. Nevertheless, the dust temperature that 
we derive is probably a good estimate of the kinetic tem- 
perature of the HCO^/SiO molecular clump surrounding 
the HC and UC Hii region (see Sect. below). 



In Table ||, for each source and at each frequency, 
we report the peak position (a^'gso and <5i95q), the peak 
^^pcak^ and integrated fluxes (F^) , the observed (Oscr) 
and deconvolved (Q-ia) angular diameter, and the corre- 
sponding linear radius (rao-) of the observed continuum 
emission (dust plus free- free. Figs. |] and || top panels). 
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Table 3. Derived parameters of the dust continuum emission from the cores 
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Fig. 4. G29.96— 0.02 source integrated spectra of the 
HCO+(1-0) and SiO(2-l) v=0 hues. The vertical Unes 
indicate the velocity ranges used for the maps shown in 
Fig. H (see text for details) . 

In Table ^ the same quantities are reported for the esti- 
mated dust emission (Figs. |^ and || bottom panels). The 
observed angular diameters of the sources have been com- 
puted as 630- = 2^fApK, where A is the area encompass- 
ing the observed source emission inside the 3 cr contour 
level. The values of 83^ have then been deconvolved us- 
ing the expression ^3^ = ^63^ — 83, where 8b is the 
synthesized HPBW. The linear radii have been calculated 
from 03o- using the adopted distances of 6 and 7.9 kpc, 
for G 29.96 -0.02 and G31.4 1+0.3 1 respectively (Pratap 
et al. 119991; Churchwell et al 



1990) 



3.2. Line 

In Figs. ^ and || we show the source integrated spec- 
tra for both HCO+(1-0) and SiO(2-l) v=0 transitions. 
The HC0"'"(1— 0) spectra show absorption at a velocity 
very close to that of the ambient cloud in the case of 
G31.41-h0.31, and redshifted in G29.96-0.02. No absorp- 
tion is revealed in the SiO(2— 1) v=0 spectra. From the 
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Fig. 5. As Fig. |, for G31.41+0.31. 

spectra, we selected the range of velocity channels high- 
lighted by the vertical bars in Figs. ^ and ^, in order to 
produce maps of the HCG+ absorption (dotted bars) and 
of the line wings (dashed for the red wings and dash-dotted 
bars for the blue wings). We produced also maps averaged 
over all the line emission. 

The contour maps of the averaged HCG+(1— 0) and 
SiO(2— 1) v=0 emission, superimposed on the correspond- 
ing continuum in grey scale, are displayed in Fig. ^ for 
G29.96-0.02 and Fig. for G31.4H-0.31. The molecular 
emission in our HCO"'" and SiG maps of G29.96— 0.02 is 
elongated toward the position of the CH3OH clump found 
by Pratap et al. ( |l999| ). Particularly, the SiO map shows a 
second peak close to the estimated po sition of the CII3OH 
outflow centre found by Pratap et al. ( 1999 ). As expected, 
in both sources the region of the HGO"*" absorption is co- 
incident with the position of the 88 GHz continuum peak, 
as shown in Fig. this proves that it is a true absorption 
feature, not a self-absorption effect. 

Towards G29.96— 0.02, the peak of the molecular emis- 
sion is coincident with the core seen in the dust continuum 



and NH3 emission by Cesaroni et al. ( 1994a ). We note 
that in this source most of the 88 GHz continuum is due 
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Fig. 6. Maps toward G29.96-0.02 in the HCO+(1-0) 
and SiO(2— 1) v=0 emission averaged in velocity over all 
the line (top) and over the blue (solid contours) and red 
(dashed contours) wings (bottom), superimposed to the 
88 GHz continuum map (grey scale). The ranges of ve- 
locity used in the maps are: for the HCO"'" 84-111 km/s 
for all the line, 84-89 km/s for the blue wing and 106- 
111 km/s for the red wing; for the SiO 85-111 km/s for 
all the line, 85-92 km/s for the blue wing and 104-111 
km/s for the red wing. The contour levels for the HCO+ 
line are: -0.10, 0.10, 0.15, 0.20, 0.25, 0.3, 0.35, 0.4, 0.43 
Jy/beam; -0.09, 0.09, 0.12, 0.15, 0.17 Jy/beam (blue) 
and -0.09, 0.09, 0.12, 0.15 Jy/beam (red). The contour 
levels for the SiO are: -0.05, 0.05, 0.10, 0.15, 0.20, 0.25 
Jy/beam; -0.09, -0.07, 0.07, 0.09, 0.12, 0.15, 0.18, 0.20 
Jy/beam (blue) and -0.09, -0.07, 0.07, 0.09, 0.12, 0.145 
Jy/beam (red). The star marks the position of the UC Hii 
region ionizing star, the circle is the 88 GHz dust peak, 
the filled diamond the 24 GHz continuum peak, the tri- 
angle the NH3(4,4) HC peak, and the cross marks the 
pea k pos ition of the CH3OH clump detected by Pratap et 
al. (119991) . 



to the free-free emission from the UC H 11 region, which is 
slightly displaced from the dust core (see Sect. 



2.1). The 



lack of HCO+ emission towards the UC H 11 peak is due to 
the absorption (see Fig. 0). The red and blue-shifted wing 
emissions of both HCO~ and SiO appear only slightly 
separated along a E-W axis, parallel to that defined by 
th e velo city gradient of the ammonia core (Cesaroni et 
al. 119981) . 

Towards 031.41+0.31, the peak of the HCO+ and SiO 
clump is shifted southwards with respect to the position 
of the core traced by the dust emission and high exci- 
tation NH3 an d CH3 CN transitions (Cesaroni et al. 1998 ; 
Cesaroni et al. 1994b). In the case of HCO+, this is proba- 




HCO*(1-0) 
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Fig. 7. Maps toward G31.4H-0.31 in the HCO+(1-0) 
and SiO(2— 1) v=0 emission averaged in velocity over all 
the line (top) and over the blue (solid contours) and red 
(dashed contours) wings (bottom), superimposed to the 
88 GHz continuum map (grey scale) . The ranges of veloc- 
ity used in the maps are: for the HCO"'" 84-111 km/s for 
all the line, 84-89 km/s for the blue wing and 106-111 
km/s for the red wing; for the SiO 83-111 km/s for all the 
line, 83-89 km/s for the blue wing and 106-111 km/s for 
the red wing. The contour levels for the HCO"*" line are: 
-0.05, 0.05, 0.10, 0.20, 0.25, 0.27 Jy/beam; -0.08, 0.08, 
0.10, 0.12, 0.14, 0.16 Jy/beam (blue) and -0.08, 0.08, 0.10, 
0.12 Jy/beam (red). The contour levels for the SiO are: 
-0.05, 0.05, 0.08, 0.12, 0.16, 0.20, 0.24 Jy/beam; -0.09, 
0.09, 0.12, 0.15, 0.18 Jy/beam (blue) and -0.09, 0.09, 0.12, 
0.14 Jy/beam (red). The square marks the position of the 
110 GHz continuum peak, the circle is the 88 GHz dust 
peak, the filled diamond the 24 GHz continuum peak and 
the triangle the NH3(4,4) peak. 



bly due to the fact that the core is coincident with the con- 
tinuum peak, where the line is in absorption (see Fig. ||). 
In this case, the red and blue wings are clearly separated 
along a NE-SW axis, parallel to the velocity gradient seen 
in NH3 and CH3CN. The SiO red wing overlaps with the 
position of the dust continuum emission. 

In Table ^ we report the parameters of the HCO"*" line 
emission. For each source we report the position, velocity, 
and synthesised beam brightness temperature (Tsb) for 
the emission and absorption peaks. In the same table we 
also report the observed (63^) and deconvolved (Osa) an- 
gular diameter of the HCO"*" molecular clu mp, c omputed 
from the 3ct contour, as explained in Sect. 3J, and the 
deconvolved radius (rsa)- In Table || we report the param- 
eters for the SiO transition: position, velocity, and synthe- 
sised beam brightness temperature (?sb) for the emission 
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Fig. 8. Contour plots of the HC0+(1— 0) line averaged 
over the absorption velocity range, 99-102 km/s for 
G29.96-0.02 and 95-99 km/s for G31.41+0.31, marked by 
the dotted lines in Figs. ^ and ||, overlaid on the 88 GHz 
continuum emission. Negative contours are dashed. The 
contour levels are -0.5, -0.4, -0.3, -0.2, -0.1, 0.1, 0.15, 
0.2, 0.25 Jy/beam for G29.96-0.02 and -0.42, -0.4, 
-0.35, -0.3, -0.254, -0.2, -0.15, -0.1 Jy/beam for 
G31.41+0.31. Symbols are as in Figs. |and@. 



peak. We also give the values of the integrated line emis- 
sion (/ TdV) and the peak velocity (I^sr) from gaussian 
fits to the spectra obtained by averaging in each spectral 
channel the line emission within the contour correspond- 
ing to 50% of the peak in Figs. ^ and ^. The observed 
and deconvolved angular size corresponding to these con- 
tours and the deconvolved linear radius are also reported 
in Table |. 



4. Discussion 

4.1. Mass estimates 

Assuming optically thin dust thermal emission at 88 GHz, 
we calculated the mass of the millimeter continuum 
sources using: Md = -;r^^^TrS^ where F^, is the observed 



K,,B„{T) ' 

integrated flux density on dust emission map (i.e. after 
subtracting the computed free- free emission), d the dis- 
tance from the Sun, B^{T) the Planck function at the 



assumed dust temperature T, = K230 GHz I _ 
the dust opacity per gram with K230 ghz = 0.005 cm^ g~^, 
assuming a gas-to-dust ratio of 100 by mass (Preibisch et 
al. 1993), and (3 is taken from Table ^ 

The gas kinetic temperatur es of th e two HCs have been 
estimated by Cesaroni et al. ( 1994a ) using the NH3(4,4) 
line; they found 83 K for G29.96-0.02 and 120 K for 
G31. 41-1-0. 31. This latter value is consistent with the es- 
timate of 123 K base d on the CH3CN observations by 
Cesaroni et al. ( 1994b| ). We used the temperature derived 
from the hot molecular gas, rather than the lower temper- 
ature derived from the global SED fit, because the latter is 
an estimate of the extended envelope temperature. In fact, 
the temperature is determined by the grey-body spectrum 
peak, i.e. the IRAS flux values, wich are sensitive to the 
extended dust emission, while the continuum emission de- 
tected in our 88 GHz maps is dominated by the hot dust 
within the HC. With these temperatures, we estimated 
the mass of the hot core responsible for the 88 GHz dust 
continuum emission: 2900 Mg for G29.96-0.02 and 1500 
Mq for G31.4H-0.31. The latter is consistent with the es- 
timate obtained from the 233 GHz flux for the same source 
(1200 Mq). While, for G29.96-0.02 due to the problems 
discussed in Sect. 2.1, the 233 GHz map cannot be used 



to derive a mass estimate. 

The total mass of the molecular clumps surrounding 
the UC Hii and HC can be estimated from the HCO+ 
and SiO transitions. Here we will attempt to evaluate the 
masses using the virial and LTE assumptions. Given the 
velocity structures and complex line shapes observed in 
both lines, these estimates are affected by large uncer- 
tainties. 

Assuming that the HC0+(1— 0) is optically thick, the 
peak TsB (Table 0) can be used as an estimate of the 
gas kinetic temperature. Another estimate of the temper- 
ature of the molecular clump as a whole can b e ob tained 
from the continuum SED fits presented in Sect. 3.1. These 
two estimates are consistent for G29.96— 0.02, while for 
G31. 41-^0. 31 the HCO+ estimate is much lower. This dis- 
crepancy probably arises from the fact that in this source 
the absorption occurs towards the molecular clump peak, 
while in G29.96— 0.02 the molecular clump, centered on 
the HC, is shifted from the position of the UC Hii where 
the absorption is centered. We decided to adopt an aver- 
age value of 30 K for both clumps. We assumed an SiO 
abundance of 7x10^^^ with respect to H2, following the 
average value found by Codella et al. ( |1999D for the quies- 
cent SiO component in a sample of young stellar objects. 
This value is also in agreement with the abundance de- 
rived by McMullin et al. (200C ) for the quiescent gas in the 
Serpens core. Using the integrated line emission given in 
Table |, we obtained a mass of 1600 M© for G29.96-0.02 
and 3200 M© for G31. 41-^0. 31. The main uncertainties 
in these estimates are the assumed SiO abundance and 
LTE-cquilibrium condition; in fact the SiO molecule abun- 
dance is well known to be enhanced by a factor 1 0-100 
in shocked regions (Bachiller & Perez Gutierrez 1997 ; 
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source 




01950 


5l950 






63. 


6'ia 












(km s^-"") 


(K) 


(") 


(") 


(pc) 


G29.96-0.02 


emission 


18:43:26.87 


-2:42:37.8 


95.5 


25.7 


15.3 


14.6 


0.21 




absorption 


18:43:27.15 


-2:42:37.1 


98.8 


-6.0 








G31.41+0.31 


emission 


18:44:59.15 


-1:16:11.5 


92.1 


17.2 


14.4 


13.7 


0.26 




absorption 


18:44:59.15 


-1:16:07.3 


97.2 


-5.1 









Table 5. Parameters of the SiO(2— 1) v=0 line emission 



source 


01950 


^1950 






Vlsr 


AVfwhm 


©50% 


^50% 


^50% 








(K) 


(K km s^^) 


(km s"^) 


(km s^^) 


(") 


(") 


(pc) 


G29.96-0.02 


18:43:26.91 


-2:42:39.2 


8.1 


39.3 


97.9 


10.5 


7.7 


6.1 


0.09 


G31.41+0.31 


18:44:58.96 


-1:16:10.8 


6.6 


39.8 


95.7 


15.3 


8.6 


7.2 


0.14 



Codella et al. 19991 ). Such abundance enhancement would 
lead to unrealistically high mass estimates, unless the SiO 
is confined to an outflow. 

To derive the virial mass from the SiO(2— 1) v=0 tran- 
sition we u sed the FWHM calculated from the gaussian fit 
(Sect. U, Table |). We estimated a mass of 2500 Mq for 
G29.96-0.02 and 8000 Mq for G31.41+0.31. Owing to the 
fact that the SiO emission might arise from outflowing gas, 
which would increase the line FWHM, these values should 
be considered as upper limits. In the case of HC0+(1— 0), 
the presence of the absorption feature makes very diffi- 
cult to determine the intrinsic FWHM of the line. We fit- 
ted the line profile with two gaussian components, one in 
emission and one in absorption. The resulting linewidths 
of the emission component are 6.1 and 7.2 kms^^ for 
G29.96-0.02 and G31.41+0.31, respectively. The result- 
ing virial mass estimates, using the sizes from Table ^ are 
2000 and 3400 A/© for G29.96-0.02 and G31.41-f 0.31, re- 
spectively. However, the fit is strongly biased by the line 
wings and thus the resulting virial mass is most proba- 
bly overestimated. This problem is particularly evident in 
G31.41-f0.31, where these features are detected in both 
SiO and HCO"'', and the virial mass estimates are less re- 
liable and a few times higher than the masses derived from 
the dust emission or using the LTE approximation. 

In conclusion, considering all the uncertainties in each 
of the methods discussed above, our estimates of the 
masses of the molecular clumps in both sources are con- 
sistent with 2000 Mq with an error of ~5G%. 



4.2. Gas kinematics 

We investigated the presence of a velocity gradient along 
the axes parallel and perpendicular to the direction where 
such a grad ient was found in the NH3(4,4) (Cesaroni 
et al. 1998 ) and CH3 CN(6 -5) transitions (Cesaroni et 
al. |l994bt Olmi et al. |l996|) . For both the HC0+(1-G) 



and SiO(2— 1) v=0 lines we averaged the emission along 
the axis from the UC H ii region peak to the core, and 
perpendicularly to it. This axis is along the E-W di- 
rection for G29. 96-0.02 and along the NE-SW direc- 
tion for 031.41+0.31. In Figs. | and |l| the position- 
velocity diagrams are shown, respectively for G29.96— 0.02 
and G31. 41-1-0. 31: we plotted, in the top panics for 
HCO+(1-0) and in the bottom for SiO(2-l) v=0, the 
offset along the axes previously described. In SiO a 
velocity-position trend is clearly seen along the "paral- 
lel" axes. The Vlsr varies from 89 to 107 km s^^ over 
-3" toward G29.96-0.02 and from 86 to 106 km s'^ 
over ~4" toward G31.41+0.31. The corresponding veloc- 
ity gradients are approximately 6 km s~^ arcsec"^ toward 
029. 96-0. 02 and 5 km s"^ arcsec"! toward G31.41-h0.31, 
consistent with the velocity gradient of approximately 
2 — 3 km s^^ arcsec"^ over 3" found along the same 



directions in ammonia by Cesaroni et al. (1995). A ve- 
locity gradient of the same order is seen in HCO+ to- 
ward 031. 41-1-0. 31, while it is not evident in the case of 
029.96— 0.02, most likely because of the presence of ab- 
sorption. In the following we propose an interpretation of 
these findings 



4.2.1. G29.96-0.02 

In the case of the 029.96— 0.02 molecular clump, as shown 
in the position- velocity diagrams (Fig. ^) , the center of the 
H0O+ clump is coincident in both position and velocity 
with the NH3 HO; which in turn is located at the geomet- 
rical center of the molecular clump, as shown in Fig. 1 of 
Kurtz et al. (2000). The H00+ absorption appears red- 
shifted with respect to the molecular clump s ystem ic Vlsr 
(as tr aced b y C^^S and O^'^O; Cesaroni et al. |l99li Hofner 
et al. 2000 ), and positionally coincident with the UC Hii 
region, located east of the emission geometrical center. 
HOO"*" channel maps in Fig. ^ also show blueshifted gas 
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Fig. 10. Channel maps of HC0+(1— 0) line towards G29.96— 0.02, overlaid on the 88 GHz continuum emission. Contour 
levels are from —0.55 to —0.15 and from 0.15 to 2.45 by 0.1 Jy/beam. Negative contours are dashed. Symbols are as 
in Fig. ^. 



emission toward the same direction, suggesting that the 
gas behind the H ii region is moving towards the observer. 
This implies that the absorbing gas is moving along the 
line of sight towards the UC H ii region. Both pv-diagrams 
also show that the maximum line-of-sight gas velocities are 
in the direction of the core position. 

All these features can be easily interpreted if the UC 
H II is embedded within an infalling, spherically symmet- 
ric HCO"'" clump centered on the HC position. As shown 
in the sketch in Fig. |l^, if the clump is contracting toward 
the NH3 HC, toward the UC H 11 we will observe redshifted 
absorption from the gas in front of it and blueshifted emis- 
sion from the gas behind it. At the position of the HC, 
there is emission at the highest line of sight velocities for 
both the redshifted and blueshifted gas. 

The less abundant, optically thin SiO molecule traces 
a more internal region of the molecular clump. For this 
molecule, the pv-diagrams show a peak coincident with 
the HC, within the uncertainties, and a velocity structure 
consistent with that derived from the NH3(4,4) observa- 
tions of Cesaroni et al. ( |1998D . To explain this feature, 
we can hypothesize that SiO could trace the outer regions 
of the rotating NH3(4,4) disk or the inner regions of the 
infalling molecular clump (see Fig. |l^). In the first sce- 
nario, the velocity gradient would be due to the velocity 
field in the outer regions of the disk. In this view the SiO 
abundance could be enhanced at the interface between 
the infalling molecular gas and the surface of the disk. 
Altenatively, the SiO emission could trace the innermost 



and denser regions of the infalling envelope: the observed 
velocity gradient could then be produced by the asym- 
metric erosion by the UC H 11. If the ionised gas is sligthly 
off-center and further away from the observer with respect 
to the center of the infalling clump (as sketched in Fig. pi] ), 
blue-shifted SiO gas is eroded from one of the sides of the 
clump giving the illusion of a velocity gradient accross the 
core. It should be noted, however, that the SiO is generally 
thought to trace outflowing material, as shown by several 



authors (Bachiller & Perez Gutierrez 1997; Cesaroni et 
al. 1999a; Acord et al. 1997), therefore one cannot ex- 



clude a priori the possibility that the velocity structure 
observed at the HC position (Fig. |6[ bottom panels) may 
be due to a compact jet-outflow component unresolved by 
our observations. 

At the p osition of the CH3OH core detected by Pratap 
et al. ( 1999 ) we detect emission in both SiO and HCO+. 
We do not detect any clear signature of outflowing gas at 
that location, and favour the interpretation of a second, 
less massive molecular clump. Future higher resolution ob- 
servations may allow to clarify the matter definitively. 



4.2.2. G31.41+0.31 

Compared to G29.96-0.02, the structure of G31.41-h0.31 
is more difficult to interpret. In this case, the bulk of the 
free- free emission from the UC Hii region is along the 
same l ine of s i ght o f the NH3(4,4) HC detected by Cesaroni 
et al. ( 1994a 1998 ). Due to the absorption at the position 
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Fig. 13. As fig. for G31.41+0.31. Tlie contour levels are from -0.45 to -0.15 and from 0.15 to 1.65 by 0.1 Jy/beam. 
Symbols are as in Fig. ^ 



of the HC, the observed molecular line emission peak, both 
in HCO+ and SiO, is south of the HC (Figs. and |). 

The pv-diagrams (Fig. |l^) show at high velocity, V > 
8 km/s from the HC Vlsr, in both transitions, a gradi- 
ent consistent with that observed in NH3(4,4) (Cesaroni 
et al. 199S). At lower velocity, however, the situation in 



HCO"*" is more complicated, with emission and absorption 
occurring toward the same line of sight at both blueshifted 
and redshifted velocities, with respect to the Vlsr of the 
HC and of the m olecul ar clump as see n in C '^^S and C^^O 
(Cesaroni et al. 1991 Hofner et al. 2000 ). The channel 
maps (Fig. |l3|) clearly show the superposition of these ve- 
locity features, and show an extension of the emission at 
blue velocities to the south of the HC position. In SiO, 
the emission peak lies south-west of the NH3 HC and is 
blueshifted with respect to the latter (Figs. and [l^) . The 
red emission (Fig. |^) peaks toward the HC, whereas the 
line bulk and the blue emission peak south-west of the 
HC, the velocity gradient agreeing with that seen in NH3. 

It seems that there is a north-south gradient from red 
to blueshifted emission in HCO+, possibly tracing an out- 
flow or the outer regions of the NH3 disk, and a NE- 
SW gradient in SiO, compatible with the disk hypothe- 
sis. Thus, while the observed velocity field could also be 
interpreted as a combination of infall, outflow and rota- 
tion, as it m ay happen in similar regions (see e.g. Wilner 
et al. 1996 ), it is clear that we lack the spatial resolution 
required to properly separate the various components. In 
fact, at present we cannot exclude the possible alternative 
explication that we are actually observing the superposi- 



tion of two molecular cores, one located at the NH3 HC 
position and the second south of it, where the HCO+ and 
SiO emission peak. In this view, the velocity structure we 
observe would be partially produced by the superposition 
of the velocities of the two cores. However, we note that 
unlike G29.96-0.02, in G31. 41-^0.31 the absorption dip 
is seen at the centre of the HCO^ line: this seems to be 
more consistent with the disk hypothesis, as collapse and 
outflow should determine a shift in velocity of such an 
absorption dip. In conclusion, we would tend to favour 
the disk hypothesis, although possibly in a more complex 
scenario involving outflow and/or infall. 

5. Conclusions 

Our interferometric observations show that the HCs and 
UC Hii in G29.96-0.02 and G31.41-H0.31 are embed- 
ded within massive gaseous clumps with a complex ve- 
locity structure. The masses of the clumps as derived 
from the millimeter continuum emission and from virial 
and LTE estimates based on the molecular line emis- 
sion are in the range 1000-3000 Mq for both sources 
(with G31. 41-1-0. 31 being sligthly more massive than 
G29.96-0.02). In G29.96-0.02 the HC is located close 
to the geometrical center of the clump and there is com- 
pelling evidence for a global motion of the molecular ma- 
terial onto the central object. In G31. 41-1-0. 31 our spatial 
resolution is not high enough to resolve the various com- 
ponents and it is not easy to disentangle the observed 
velocity structure. Most likely, in this source several kine- 
matical components (outflow, rotation and possibly infall) 
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Fig. 9. Contour plots for G29.96-0.02 as a functfon of 
Vlsr and offset from the phase center afong the E-W axis 
from the Hii region and the core peak (left) and perpen- 
dicularly to it (right). The maps have been obtained by 
averaging the HCO+(1-0) (top) and SiO(2-l) v=0 (bot- 
tom) emission in G29.96-0.02 over 8'.'4 along the E-W 
axis and over 7" perpendicularly to it. The contour levels 
are: for HCO+, in Aa -0.4, -0.35, -0.25, -0.15, -0.09, 
0.09, 0.15, 0.25, 0.50, 1.0, 1.5, 2.0 Jy/beam, in A6 -0.09, 
0.09, 0.15, 0.25, 0.50, 1.0, 1.5 Jy/beam; for SiO, in Aa 
-0.15, -0.09, 0.09, 0.15, 0.20, 0.30, 0.40, 0.50 Jy/beam, 
in A6 -0.09, 0.09, 0.15, 0.20, 0.40, 0.49 Jy/beam. The 
dot-dashed lines mark the position and the velocity of the 
core detected in NH3(4,4) by Cesaroni et al. ( |1998D . 



are merged together by the spatial resolution of our ob- 
servations. 

In both regions our observations are consistent 
with earlier claims of the presence of rotating disks 
around young accreting massive (proto-)stars (Cesaroni et 
al. 1998). However, our spatial resolution does not allow 
us to definitely prove such claims and alternative models 
to interpret our observations are still viable. To definitely 
solve this dilemma and to probe the detailed structure of 
G31.41-f0.31, higher resolution observations are needed. 
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the gas seen in absorption toward the UC Hii. The ob- 
server is on the left side. 



Afflerbach A., Churchwell E., Hofner P., Kurtz S., 1994, ApJ, 
437, 697 

Bachiller R., Perez Gutierrez M., 1997, ApJ, 487, L93 

BaU R., Meixner M.M., Keto E., Arens J.F., Garrett Jernigan 

J., 1996, AJ, 112, 1645 
Cesaroni R., Felli M., Jenness T., Neri R., Olmi L., Robberto 

M., Testi L., Walmsley CM., 1999a, A&A, 345, 949 
Cesaroni R., Felh M., Walmsley CM., 1999b, A&AS, 136, 333 
Cesaroni R., Felh M., Testi L., Walmsley CM., Olmi L., 1997, 

A&A, 325, 725 

Cesaroni R., Hofner P., Walmsley CM., Churchwell E., 1998, 
A&A, 331, 709 

Cesaroni R., Walmsley CM., Kompe C, Churchwell E., 1991, 
A&A, 252, 278 

Cesaroni R., Churchwell E., Hofner P., Walmsley CM., Kurtz 

S., 1994, A&A, 288, 903 
Cesaroni R., Olmi L., Walmsley CM., Churchwell E., Hofner 

P., 1994, ApJ, 435, L137 
Chini R., Kreysa E., Mezger P.G., Gemiind H.-P., 1986, A&A, 

154, L8 

Churchwell E., Walmsley CM., Cesaroni R., 1990, A&AS, 83, 
119 



C. Maxia et al.; The hot cores in G29. 96-0.02 and G31.41+0.31 



13 



G31 .41 +0.31 

A|| (orcsec) A| (arcsec) 

10 -10 -10 10 




10 -10 -10 10 



A|| (arcsec) (crcsec) 

Fig. 12. Contour plots for G31. 41-1-0. 31 as a functfon of 
Vlsr and offset from the phase center along the SW-NE 
axis from the H ii region and the core peak (left) and per- 
pendicularly to it (right). The maps have been obtained 
by averaging the HCO+(1-0) (top) and SiO(2-l) v=0 
(bottom) emission in G31. 41-1-0. 31 over 10'.'5 along the 
SW-NE axis and perpendicularly to it. The contour levels 
are: for HCO+, in A|| -0.30, -0.25, -0.20, -0.15, -0.10, 
-0.06, 0.06, 0.10, 0.25, 0.50, 0.75, 1.0, 1.12 Jy/beam, in 
A_L -0.3, -0.25, -0.20, -0.15, -0.10, -0.06, 0.06, 0.10, 
0.25, 0.50, 0.75, 1.0 Jy/beam; for SiO, in A,, -0.1, -0.06, 
0.06, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.38 Jy/beam, in 
A_L -0.09, -0.06, 0.06, 0.10, 0.15, 0.20, 0.25, 0.30, 0.32 
Jy/beam. The dot-dashed lines mark the position and the 
velocity of the core detected in NH3(4,4) by Cesaroni et 
al. ( 119981) . 

Codella C, Bachiller R., Reipurth B., 1999, A&A, 343, 585 
Condon J.J., Cotton W.D. et al., 1998, AJ, 115, 1693 
Fey A.L., Gaume R.A., Claussen M.J., Vrba F.J., 1995, ApJ, 
453, 308 

Gaume R.A., & Mutel R.L., 1987, ApJS, 65, 193 

Herter T. et al., 1981, ApJ, 250, 186 

Hofner P., ChurchweU E., 1996a, A&AS, 120, 283 

Hofner P., Kurtz S., ChurchweU E., Walmsley CM., Cesaroni 

R., 1996b, ApJ, 460, 359 
Hofner P. Kurtz S., ChurchweU E., Watson A.M., 1996c, 

BAAS, 188, 17 

Hofner P., Wyrowski F., Walmsley CM., ChurchweU E., 2000, 

ApJ, 536, 393 
Hunter T.R., Phillips T.G., 1997, ApJ, 478, 283 
Kurtz S., Cesaroni R., ChurchweU E., Hofner P., Walmsley 

CM., 2000, in Protostars & Planets IV, ed. V. Mannings, 

A. P. Boss, & S.S. Russell (Tucson: Univ. Arizona Press), 

p. 299 

McMulhn J.P., Mundy L.G., Blake G.A., Wilking B.A., 
Mangun J.G., Latter W.B., 2000, ApJ, 536, 845 

Molinari S., Brand J., Cesaroni R., Palla F., 1996, A&A, 308, 
573 



Molinari S., Testi L., Brand J., Cesaroni R., Palla F., 1998, 
ApJ, 505, L39 

Mooney T., Sievers A., Mezger P.G., Solomon P.M., Kreysa 

E., Haslam C.G.T., Lemke R., 1995, A&A, 299, 869 
Olmi L., Cesaroni R., Walmsley CM., 1996, A&A, 315, 565 
Palla F. & Stahler S.W., 1993, ApJ, 418, 414 
Pratap P., Menten K.M., Snyder L.E., 1994, ApJ, 430, 129 
Pratap P., Megeath S.T., Bergin E.A., 1999, ApJ, 517, 799 
Preibisch Th., Osseukopf V., Yorke H.W., Henning Th., 1993, 
A&A, 279, 677 

Scoville N.Z., Carlstrom J.E., Chandler C.J., Phillips J. A., 
Scott S.L., Tilanus R.P.J. , Wang Z., 1993, PASP, 105, 1482 
Testi L., Felli M., Persi P., Roth M., 1998, A&A, 329, 233 
Testi L., Hofner P., Kurtz S., Rupen M., 2000, A&A, 359, L5 
Tofani C, Felli M., Taylor G.B., Hunter T.R., 1995, A&AS, 
112, 299 

Walsh A.J., Burton M.G., Hyland A. P., Robinson G., 1998, 

MNRAS, 301, 640 
Watson A.M., Hanson M.M., 1997, ApJ, 490, L165 
Watson A.M., CoU A.L., Shepherd D.S., Hofner P., ChurchweU 

E., 1997, ApJ, 487, 818 
Wilner D.J., Ho P.T.P., Zhang Q., 1996, ApJ, 462, 339 
Wink J.E., Altenhoff W.J., Mezger P.G., 1982, A&A, 108, 227 
Wood D.O.S., ChurchweU E., 1989a, ApJ, 340, 265 
Wood D.O.S., ChurchweU E., 1989b, ApJS, 69, 831 
Wood D.O.S., ChurchweU E., Salter C.J., 1988, ApJ, 325, 694 



